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Abstract 

In this paper we study the one dimensional dynamical Casimir effect. We con- 
sider a one dimensional cavity formed by two mirrors, one of which performs 
an oscillatory motion with a frequency resonant with the cavity. The naive 
solution, perturbative in powers of the amplitude, contains secular terms. 
Therefore it is valid only in the short time limit. Using a renormalization 
group technique to resum these terms, we obtain an improved analytical so- 
lution which is valid for longer times. We discuss the generation of peaks in 
the density energy profile and show that the total energy inside the cavity 
increases exponentially. 

I. INTRODUCTION 

The problem of quantum fluctuations inside cavities has attracted the attention since a 
long time ago [HI]. A way to study the structure of the vacuum is to distort it by changing 
the configuration of the cavity in time PJ^. The simplest configuration is that of a one 
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dimensional cavity formed by two perfectly reflecting mirrors, one of which is fixed and the 
other one is allowed to move in a predetermined way, or rather, its motion is determined by 
the backreaction of the electromagnetic field [Q . There are in fact only a few predetermined 
motions which allow an exact resolution of the problem. In |^ a special motion for the 
mirror, which has an exact solution, has been considered, and it has been shown that the 
Casimir force may be resonantly enhanced. In |P] a geometrical method for solving the 
problem for arbitrary wall motions has been developed, and basically the same structure for 
the electromagnetic field within the cavity has been found. Of special interest are the cases 
where the moving mirror oscillates with one of the eigenfrequencies of the unperturbed cavity 
1^-^. A naive approach is to make perturbations in the amplitude of oscillation. However, 
this perturbative treatment has only a very limited range of validity: the appearance of 
secular terms proportional to the time imply that after a short period the approximation 
breaks down. In this paper we apply a method inspired in the Renormalization Group (RG) 
to treat these singular perturbations. The method has a wide range of application ||10|| , 
especially to ordinary differential equation problems involving multiple scales, boundary 
layers, asymptotic matching and WKB analysis. The main advantage of the RG method 
is to provide a simple and unified calculational method for all problems of this sort. In 
our present case, the application of the method to the dynamical Casimir effect permits 
to get a solution for the structure of the electromagnetic field within the cavities that is 
valid for a period of time longer than that of the perturbative case. With the solution at 
hand, we study local properties of the field such as the energy density. In agreement with 
other authors 0, we show that the resonant moving wall induces an exponential growth of 
the total energy, and that peaks form inside the cavity, which travel at the speed of light 
bouncing against the walls. 

We consider a one dimensional cavity formed by two perfectly reflecting mirrors. One 
of them is fixed at x = 0, while the other one performs an oscillatory motion L{t) = 
Lq[1 + esin(^)], with e -C 1 and g G N, i.e., the moving mirror oscillates with a frequency 
equal to one of the eigenfrequencies of the cavity. We shall assume that the oscillations 
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begin at t = 0, and that the mirror is at rest for t < 0. Note that we shall not treat the 
moving mirror as a degree of freedom (either classical or quantum), but just as a given time- 
dependent boundary for the electromagnetic field inside the cavity. The vector potential 
A{x,t) satisfies the one dimensional field equation DA = and the boundary conditions 
A{x = 0,t) = A{L{t),t) = for all times. Therefore one can express the field inside the 
cavity as 



oo 

A{x,t) = ^^aktpk{x,t) + altjjl{x,t) , (1) 



k=l 



where the mode functions ipk{x, t) are positive frequency modes for t < 0, and and are 
time-independent annihilation and creation operators, respectively. 
If one writes the modes in terms of a function R(t) as [@] 

the boundary conditions are met provided that 

R{t + Lit)) - R{t - Lit)) = 2. (3) 

The complete solution to the problem involves finding a solution R{t) in terms of the pre- 
scribed motion L{t). The modes are positive frequency modes for t < if R{t) = t/L^ for 
—Lq < t < Lq, which is indeed a solution to Eq.(|^) for t < 0. Note that the boundary 
condition for R{t) involves its values over the whole range of time —Lq <t< Lq. 

In what follows, we will describe a method to find an analytical approximation to the 
solution of Eq.(|]). We shall first obtain a perturbative solution by expanding in powers 
of the amplitude e (Section 3). As this perturbative solution will contain secular terms, it 
will be valid only for short times, i.e. e < 1. Using RG techniques, we will be able to 
perform a resummation of the secular terms and obtain an analytical approximation valid 
for a longer period of time e^-^ < 1 (Section 4). In Section 5 we will use this solution to 
describe the evolution of the mean value of the energy density inside the cavity. In order 
to get acquainted with the renormalization group method, in Section 2 we will describe a 



simple example where the resummation of the secular terms is performed for a particular 
ordinary differential equation. The reader already familiarized with the work of Ref. 
can skip the next Section. 



II. A SIMPLE EXAMPLE 



Let us consider the Rayleigh equation 




) 
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dy 
dt 




(4) 



where e is a small number. This is an interesting oscillator because it can be shown that, for 
any initial condition and any positive e, the exact solution becomes periodic at long times 



The Rayleigh equation can be solved perturbatively using an expansion in powers of e, 
that is, y = Vo + ^Vi + C^(e^)- Up to first order in e, the perturbative solution reads 



where Yq and Bq are constants determined by the initial conditions at arbitrary t = tQ. This 
perturbative solution does not become periodic and, therefore, it is not a good approximation 
for long times. Indeed, due to the presence of the secular term, the naive perturbative 
solution is valid only for times close to the initial time to, and breaks down for e (t — to) > 1- 
This is typical of systems showing parametric resonance. Usually one has a system weakly 
coupled to an external resonant force, and as one tries to make a perturbative analysis, the 
corrections possess secular terms, i.e. terms that grow linearly with time. In the following 
we shall adopt the RG method to treat singular perturbations and we shall use it to make an 
improvement to the perturbative solution, which shall be valid for longer times, e^(t— to) < 1- 
The basic idea |]TU[ is to introduce an arbitrary time r, split t — to as {t — r) + r — to, and 
absorb the terms proportional to r — to into the 'renormalized' counterparts Y{t) and 0(r) 



and therefore approaches a limit circle in phase space |TT|. 



y{t) = Fq sin(t + Go) + e 






4 



of the 'bare' parameters contained in the zeroth order solution, that is, Yq and Gq. Using 
this idea one ehminates the secular terms proportional to r — to, and the function y{t) takes 
the form 

y{t) = Fsin(t + 6) + e 1^ (^1 - (t - r) sin(t + 6) + 

^ [cos(3(t + 6)) - cos(t + e)]| + 0{e'), (6) 

where now Y and B are functions of r. Now comes the crucial point. As r does not appear 
in the original equation nor in the initial conditions, the solution Eq.(|^) should not depend 
on T. Therefore, the partial derivative with respect to r should vanish, i.e., {dy/dT)^ = 
for any t. This is the RG equation, which implies 

The solutions to these equations are 



Y{t) = Y{to) 



4 

e(r) = e(to) + 0(e'(r-to)), (8) 



where l^(to) and 0(to) are constants to be determined by the initial conditions. We still have 
the freedom to choose the arbitrary time r. The obvious choice is t = r, since in this way 
the secular term proportional to t — r in Eq. (P) disappears. Assuming the initial condition 
?/(to) = 0, ?/(to) = 2a, with a any real number, we find Y{tQ) = 2a and 0(to) = "~^o- Finally, 
the RG improved solution reads 

Y(t)^ 

y{t) = Y{t) sin(t - to) + e4-^[cos(3(t - to)) - cos(t - to)] + 0{e^), (9) 

9d 

that is valid for e^{t — to) < 1. Note that the improved solution becomes periodic and 
approaches a limit circle of radius 2 for e{t — to) ^ 1. 

It is interesting to remark the analogy with the usual RG approach in quantum field the- 
ory: to plays the role of the ultraviolet cutoff (although there are no divergences but secular 



terms here), Yq and Go are the bare couphng constants, and Y and are the renormahzed 
counterparts evaluated at the 'scale' r. The equation is 'renormalizable' because the secular 
terms can be absorbed into the bare parameters. As anticipated, the RG is a straightfor- 
ward method to obtain, from the naive perturbative solution, an improved solution which 
is valid for longer times. In this particular example, it can be shown that the RG method is 
equivalent to multiple-scale analysis [|r^ , with the additional practical advantadge that it is 
not necessary to know a priori the multiple time scales. 



III. PERTURBATIVE SOLUTION TO THE DYNAMICAL CASIMIR EFFECT 

We will now solve Eq. (|^) using a naive perturbative expansion. We expand the function 
R{t) in terms of the small amplitude e and retain first order terms only, R{t) = RQ{t)+eRi{t). 
Equating terms of the same order we get 



Ro{t + Lo) - Ro{t - Lo) = 2, 



(10) 



i?i(t + Lo) -Ri{t- Lo) = -Lo sm{^)[R',{t + Lo) + R'o{t - Lo)]. (11) 



The general solution to Eq.(|lOD is 



Ro{t) = a+ — + J2 



n>l 



. . ,mit. ^ .nnt. 
An sm(— ) + Bn cos(— ^ 
^0 J^o 



(12) 



where a. An and Bn are constants determined by the boundary condition, that is, by the 
value of R{t) for —Lo <t< Lq. Introducing this solution into Eq. ([Tl|) we obtain 



I [R,it + Lo) -R,{t- Lo)] = sin A + ^ E^(-l)" >< 

° n>l 

(g + n)TTt 



An 
Br, 



sm 



cos 



Lo 

{q + n)T{t 

L'o 



+ sm 



cos 



(g — n)7rt 
(g — n)Txt 

To 



+ 



(13) 



whose general solution reads 
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Rlit) = (-1) 



sm 



cos 



(g + n)-nt 



— cos 



n>l ^ 



sm 



(g + n)7r)f: 



+ sin 



(g — ri)7rt' 



+ 



(14) 



where g{t) is an arbitrary periodic function with period 2Lq. We see that, as in the case of 
the Rayleigh oscillator, the perturbative correction contains secular terms that grow linearly 
in time. Therefore, this approximation will be valid only for short times, that is, e < 1. 

If we assume that the boundary condition for R{t) is already satisfied by Ro{t), then 
the periodic function g{t) must be chosen in such a way that Ri{t) = for —Lq < t < Lq. 
Therefore 



gi2pLo + z) 
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sm 



nx 



Ar. 



sm 



cos 



(g + n)-nz 
(g + n)Tiz 



n>l 



+ sin 



cos 



(g — n)'KZ 
(g — n)T{z 



+ 



(15) 



where t = 2pLo + z, p = 0,1,2, and —Lq < z < Lq. Given t, the value of the integer p 
is obtained as p = |int(t/Lo) or p = |[int(t/Lo) + 1] for int(t/Lo) even or odd, respectively. 
Note that during the first period {p = 0), g{t) makes Ri{t) vanish identically. As we 
have already seen, since the mirror is at rest for t < 0, we must impose R(t) = I/Lq for 
—Lq < t < Lq. Therefore a = An = Bn = 0, and the perturbative solution reads 



t 

Tq 



t qnt 
— sm(— ) 

^0 ^0 



z ,qTiz, 
l-sin(— . 



(16) 



The naive perturbative solution to the dynamical Casimir effect has been previously dis- 
cussed in Ref. [0. However, in that work the periodic function g{t) was taken equal to zero. 
Thus, the solution obtained there does not satisfy the correct boundary condition. 



IV. RENORMALIZATION GROUP IMPROVEMENT 



We will now adapt the KG method of Section 2 in order to obtain a solution to the Eq. 
which is valid beyond the short time limit. Let us introduce the arbitrary time r and split 
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t as t — T + T. The perturbative solution can then be written as (see Eqs.(|T2D and (0)), 
R{t)=a{T)+Y, 



n>l 



nnt .nnt, 
AJt) smf -— ) + BJt) cos( -— ^ 



t - T 

BJt) 



sm 



E 

n>l 



7i(-ir < AJt 



sm 



t-T 

{q + ?2)7rt 



sm 



(g — n)7rt 



+ 



iq + n)Tit\ f iq — n)-nt 

cos I ^ — - cos ' 



+ ^7(t,r) + 0(e2), 



(17) 



J \ Lq 

where the bare parameters a, An and -B„ have been replaced by their renormalized counter- 
parts a{r),An{T) and Bnir). Here g(t,T) denotes the function g{t) of Eq. ( PTSP with the 
same replacement. Note that then g{t, r) is no more a periodic function. 

The RG equation {dR/dT)t = consists in the present case of three independent equa- 
tions 



OT Lq 

dAJr) 

= e- 



dr ' Lq 
dBjr) 7r(-l)«+i 



dr 



2Ln 



Sng + -{\n- q\A\n-g\ - {n + q)An+q} + C(e^), 
[|n - + {n + q)Bn+g] + O(e'), 



(18) 
(19) 
(20) 



where we recall that the index n is a positive integer. The solution to Eq . ([18|) is trivial: 
a(r) = ^ + K, with k a constant to be determined. If one writes An = An — A_n and 
Bn = Bn — -B-n {n > 1), whcrc the new variables satisfy 



dAni 
'dr*' 
dBm 
'dr*' 



-6mq + (m - q)A 



TT 



m—q 



(m + q)A 



m+q 



[m — q)B, 



m—q 



{m + q)B„,+g + 0{e^), 



(21) 
(22) 



then An and Bn satisfy Eqs.(|T9|) and (|20| ) respectively. Here we have introduced a new 
time T* = re7r(— l)''+^/(2Lo). Since this set of first order differential equations ensures 
the independence of the solution R{t) with r, one can set t = t, which makes the terms 
proportional to r — t in Eq. (p!7|) vanish identically. 

The initial conditions for these differential equations are dictated by the perturbative 
solution: a(0) = ^^(0) = Bm{0) = 0. This means that k = and that Bm{t) = for all t. 
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The coefficients Am are not all zero due to the presence of the inhomogeneous term ^Smq- In 
order to solve the equation corresponding to these coefficients we introduce the generating 
functional F(s,r*) = J2m^"^^"i{T*)- Using Eq.(|2T|) we see that it satisfies the following 
differential equation 



OF 



dF 



(23) 



dr* TT ds 

with boundary condition F{s, r* = 0) = 0. We make the following Ansatz for the solution 
F{s,T*) = $[e^^ 9{^)] + h,{s), where $[. . .], g{s) and h{s) are functions to be determined. 
Introducing this form of the generating functional into the differential equation, one de- 
termines the last two functions. The function $ is determined once the initial boundary 
condition is imposed. Finally the solution reads 



F{s, 



In 



nq 



(24) 



In order to get the coefficients Am we expand this solution in powers of s. In this way 



we obtain that the only non- vanishing coefficients are Am=o{'^* 



- ln(cosh gr*) and 



Am=qj = ~j tanh"'(gr*) with j G N. Note in particular that Am<o = 0, which then means 
that the original coefficients An are equal to the A^s. 

The RG- improved solution for R{t) can be obtained from Eq. (p!7|) by setting t = t. It is 
given by 

t ^ . . ,qjTrt^ 



R{t) 



(25) 



Using the explicit form of the coefficients Am we find 



J]A,,(t) sin( 



Tm In 



iqirt 



(26) 



Lq -nq 

where we have defined ^ = exp[^^^i^^j^^]. The (now non-periodic) function g{t,t) can be 
easily evaluated 



-1)'?— sm(— ^ 



^Aqj{t)qj COS 



qjTTZ^ 



2e 



i+e + (i-e)cos(^ 



(27) 



Finally, the RG-improved solution reads 



t 2 

R{t) = Imln 1 + ^ + (1 - Oe ^0 

Lo nq 



+ 



2e 



i+e + (i-e)cos(^ 



(28) 



It is worth mentioning that this solution is valid as long as e -j^ < 1, that is, the range of 
validity of the solution is longer than the perturbative one (e < 1). In figure 1 we plot 
this function for the particular case g = 4. 

A solution to Eq.(^) in the long time limit was already obtained in [§] using a different 
procedure. It coincides with our first two terms in Eq. (^). There is perfect agreement 
between both solutions at long times because it can be shown that the third term in Eq. (pH]) 
is negligible in this limit (see Appendix A). However, as we have mentioned in the previous 
Section, this term (that comes from the periodic function g{t)) is crucial for the solution to 
satisfy the correct boundary condition at short times. 



V. ENERGY DENSITY INSIDE THE CAVITY 

In order to study the local properties of the electromagnetic field inside the cavity, we 
concentrate ourselves in the energy density of the field 



{Too{x,t)) 



dA{x,t) 
di 



) + ( 



/dA{x,t) 
\ dx 



) 



(29) 



where the expectation values are taken with respect to the vacuum state. Using the well- 
known point splitting method to regularize the divergence appearing in the energy density 
0], one can obtain the following expression for the renormalized energy density, {Tqq{x, t)) = 
—f{t + x) — fit — x), where 

1 



/ 



2Att 



2\R' 2 ^ 



r\2 



(30) 



This expression involves second and third derivatives of R{t). As R'{t) is discontinuous at 
t = (2p + 1)Lq, p = 0, 1,2, . . . (see Eq. (^)), then the energy density will develop a delta 
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function singularity which will be infinitely reflected back and forth between the mirrors. 
The physical origin of this singularity is the initial discontinuity of the wall velocity. We 
will ignore this singularity in what follows. 

The structure of the electromagnetic fleld within the cavity for our solution R{t) is similar 
to that for other existing solutions in the literature. In particular, for g > 2 the energy 
density grows exponentially in the form of q travelling wave packets which become narrower 
and higher as time increases. The total energy within the plates increases exponentially at 
the expense of the energy needed to keep the plate moving. In figure 2 we show the energy 
density profile between plates for a fixed time and for the case g = 4. As time evolves, the 
peaks move back and forth bouncing against the mirrors. The height of the peaks increases 

27rqet irqet 

as e ^0 and their width decreases as e~~^, so that the total area beneath each peak, and 
hence the total energy, grows as e~ . Apart from this exponential growth, there are 'sub- 
Casimir' regions: between the peaks the energy density takes values times smaller than 
the static Casimir case, (^oo) static ~ ^2^?- '^^^ prove all these properties analytically 
by computing the energy density with the solution given in Eq. (^) and its derivatives (see 
Appendix B) . In figure 3 the energy density is shown as a function of time at the mid point 
between plates, also for the g = 4 case. 

A rather different picture appears when one considers the g = 1 case, that corresponds 
to an oscillation frequency equal to the fundamental frequency of the cavity. In this case 
the energy density at a given point oscillates in time around the static Casimir value, and 
its time average coincides with that value. 

VI. CONCLUSIONS 

In this paper we have studied the one dimensional dynamical Casimir effect of a resonant 
oscillating cavity. For this one dimensional case, the modes of the electromagnetic field can 
be expressed in terms of the solution to the so called Moore equation. We have used a 
renormalization group improvement of the naive perturbative solution and we have succeeded 
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in obtaining an analytic solution which is valid up to times t < Loe"^, thus extending the 
range of validity of the perturbative solution (t < Lge^^). We have calculated the energy 
density inside the cavity and we have shown that a non trivial structure appears, with a 
series of peaks that grow exponentially in time and move back and forth bouncing against 
the mirrors. Although this structure has already been found in a previous work [0, here we 
have presented an analytic derivation based on the renormalization group method described 
|lO| . We expect this method to be useful to analyze the more realistic situation of a 



m 



three dimensional oscillating cavity. This analysis can be performed by studying the set of 
differential equations satisfied by the modes of the electromagnetic field This topic will 
be the subject of further investigation. 
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APPENDIX A: SHORT TIME AND LONG TIME BEHAVIOUR OF R(t) 



In this Appendix we analyze the short time (e ^ 1) and long time ^ ^ 1) behaviour 
of the RG-improved function R{t) given in Eq. (pHD . Let us first split the solution as R{t) = 
Rs{t) + -Rnp(^), where 



t 2 

RJt) = Imln 

Lo nq 



iqirt 



t 2 

arctan 

Lo nq 



i?np(t)=6(-l)''f Sin(^; 

with 



2e 



-ly+'nqet 

C = exp 



(Al) 
(A2) 

(A3) 



12 



The function i?np stems from the RG-improvement of the periodic function g{t), and it 
is non-periodic. The variable z {—Lq < z < Lq) is given in terms oi t as z = t — 2pLo with 
p = 0, 1, 2, . . .. This integer p is obtained from the value that t takes as p = |int(t/Lo) or 
p = |[int(t/Lo) + 1], for int(t/Lo) even or odd, respectively. 

For the short time limit t <^ e^^Lo, these functions are 

Kp(t)^(-irf sin(^), (A4) 

which then leads to the perturbative solution given in Eg . (p^) . 

For the long time limit t ^ e^^Lg (but t < e^^Lo since this poses the upper limit for 
the validity of our RG solution), we analyze Rs and i?np separately. We want to show that 
in this limit, the latter function is negligible. This can be graphically verified, but here we 
present an analytical demonstration. The function has a first term, linear in time, and 
a second one, that for late times becomes an oscillating function. The amplitude of the 
oscillations is independent of e. Due to this second term, Rs develops a staircase form for 
long times, as shown in Figure 1. Within regions of t between odd multiples of Lq (i.e. in 
each period p), there appear q jumps, located at values of t satisfying cos(^) = =Fl, the 
upper sign corresponding to even values of q and the lower one to odd values of q. Next we 
calculate the first derivative of Rs- Since d^/dt is proportional to e^, one can differentiate 
the function Rs with respect to time, treating ^ as a constant. The first derivative is then 
R'sit) = 2^^(t), where 

m = -T- (A5) 



i + e + (i-acos(f;^) 

Using Eg. ( |A3| ) , we see that C, vanishes (diverges) exponentially for q even (odd) at long 
times. For even g, the first derivative R's develops peaks for times t„ = l^ILliLli^^ with n an 
integer. The height of these peaks grows exponentially as Between peaks, R's vanishes 
exponentially for long times. In a similar fashion, for odd q, the first derivative develops 
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peaks at t„ = In the following, we shall consider only the even case, the odd one being 
completely similar. 

Let us now analyze the function i?np- Treating once again ^ as a constant when differen- 
tiating with respect to time, this function can be expressed in terms of the first derivative 
of Rg as follows, 

R^^{t) = {-iyezsmAR',{t). (A6) 

We see that R^p is the product of a bounded factor {\z\ < Lq) times a function F(t) = 
sin(^)i?g(t) that is proportional to the first derivative of R^ and might thus be unbounded. 
We shall now show that this is not the case. Far from the position of the peaks, F is bounded 
because R'^ is. In a surrounding of t„, we express F as 

where S = {t — tn)/LQ. Firstly, we note that this function vanishes for 6 = 0, i.e. Rnp{t = 
t„) = 0. Secondly, this function has extrema equidistant from t„ located at 6± = ± — — 

and at these points F{6±) = =F Since for long times ^ — 0, we conclude that F{6±) 

is bounded by 1/Lq. Consequently, R^p is a correction of order e to the second term of 
Eq.(P). 

APPENDIX B: STRUCTURE OF THE ELECTROMAGNETIC FIELD 

In this Appendix we study briefly the structure of the electromagnetic field within the 
cavity as given by Eqs. (|29| , |30D for our solution R{t), in order to understand the form of 
the energy density profile shown in figures 2 and 3 in the long time regime > 1. As in 
Appendix A we will split R{t) as R(t) = Rg + Rnp (see Eqs.( [AT] , [A^ )). In order to analyze 



the energy-momentum tensor, we need to study the first three derivatives of the solution 
R{t). For Rs we get 

Kit) = 2^m, (Bi) 
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R':{t)= 2^1 -e)nqsm Ant), (B2) 
Kit) = 2^(1 - e)M'[il + ^') cos A + (1 - e)[l + sm\^)mt) (BS) 

m = —r ^ -J, (B4) 

Lo [i+e + (i-acos(f )_ 

where again we treat ^ as a constant. For even q, all these derivatives develop peaks at times 
tn = ^2l!d:l}Lo.^ As can be easily seen from the above equations, the height of the peaks for 
the m-th derivative of Rg is proportional to Using the same methods as in Appendix 
A, one can show (after some algebra) that near the times t„, the m-th derivative of R^p 
also has peaks whose heights are proportional to e^^*". Since e -C 1, it means that at long 
times all the derivatives of R^p are negligible with respect to those of Rg. The form of the 
energy-momentum tensor will be governed just by the first part of our solution, namely by 
Rg. 

Let us concentrate only on the contribution to the energy density which is proportional 
to R'^ (see Eq. (^O])). From the above discussion it is clear that Tqq will develop peaks which 
grow as e . Their width decreases exponentially as e so the total area of the peaks 
grows exponentially. The same holds for the total energy in the cavity. The analysis of the 
other two terms of Eq. ([30| ) lead to the same conclusion. 

The case q = 1 shows a different behaviour. Indeed, when the energy density of the 
field is computed using the derivatives of -Rs given above, there is a cancellation between 
the different contributions in Eq.(|5IlD and the final answer coincides with the static Casimir 
value. 
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FIGURES 
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FIG. 1. R{t) vs. t/Lo as given by Eg. (p8[) . The values of the parameters are q = A and e = 0.01. 
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FIG. 2. Energy density profile between plates for fixed time t/Lg 
The amplitude coefficient is e = 0.01. 



20.4 for the q = A case. 
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FIG. 3. Energy density as a function of time for the mid-point x/Lq = 0.5 between plates. The 

parameters are q = 4 and e = 0.01. 
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